Electron-phonon bound states in graphene in a perpendicular magnetic field 
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The spectrum of electron-phonon complexes in a monolayer graphene is investigated in the pres- 
ence of a perpendicular quantizing magnetic field. Despite the small electron-phonon coupling, usual 
perturbation theory is inapplicable for calculation of the scattering amplitude near the threshold 
of the optical phonon emission. Our findings beyond perturbation theory show that the true spec- 
trum near the phonon emission threshold is completely governed by new branches, corresponding to 
bound states of an electron and an optical phonon with a binding energy of the order of auio where 
a is the electron-phonon coupling and too the phonon energy. 



Introduction - Observation of the quantum Hall ef- 
fect [H [2] in graphene [3] has stimulated extensive in- 
vestigations of the quantized motion of massless Dirac 
fermions in a perpendicular magnetic field [3HH] • Partic- 
ularly, electron-phonon interaction phenomena have be- 
come a subject of active theoretical [9HI2] and experi- 
mental [I3HI6] study. The density of states of both elec- 
trons in Landau levels and optical phonons with their 
weak dispersion shows strong delta- function peaks, which 
result in sharp resonance structures in the absorption 
spectra, tunable with the magnetic field and useful for 
future device applications. 

A unique feature of the graphene single-particle spec- 
trum is that the carrier Fermi velocity vp sa 1.15 x 10 6 
m/s is independent of its energy. As a result, the appli- 
cation of a perpendicular magnetic field B quantizes the 
electronic Dirac spectrum into nonequidistant chiral Lan- 
dau levels, — fiy^nujB, which include also a character- 
istic zero-energy state [T7lfT9] . Here fi = ±1 denotes the 
electron chirality and n — 0,1,2,... the Landau quan- 
tum number, cub = V^hvp/is is the magnetic energy 
and £b the magnetic length that determines the localiza- 
tion scale of the electron cyclotron motion. This picture 
of Landau quantization in graphene is well established, 
e.g. by direct observation using scanning tunneling spec- 
troscopy [THUDS]. 

In addition to these Landau levels, recently Ref.HZOlhas 
proposed a new sequence of discrete states in graphene 
at energies e^ n ±ujQ, corresponding to the emission or ab- 
sorption of optical phonons with energy ujq- According 
to this concept, based on perturbation theory, phonon- 
assisted electron transitions between two different Lan- 
dau levels can manifest themselves as peaks in the density 
of states at frequencies £^ n ' — ± ^o- These tran- 
sitions involve optical phonons, which exist but do not 
interact with electrons. In the present paper we show 
that even for small electron-phonon coupling, a, the true 
spectrum in graphene in a perpendicular quantizing mag- 
netic field develops a fine structure in the neighborhood 
of the phonon emission threshold, e° = + Uq, 

due to the formation of new spectral branches, corre- 
sponding to bound states of an electron and of an opti- 
cal phonon. Therefore, electronic transitions should in- 



volve these bound states and should occur not at the 
phonon emission threshold energy £° n but above and be- 
low it at the characteristic scale of the binding energy, 
W gr = auj (see Fig. [TJleft)) [2T]. The bound states of 
electron and optical phonon were first proposed for bulk 
semiconductors in Ref. [22] and it is important to note 
that no resonance situation is required for the formation 
of the bound states. In graphene the electron-phonon 
binding energies are of both signs because no contin- 
uum exists in the electron+phonon system in the pres- 
ence of the magnetic field and the bound states appear 
both below and above the phonon emission "threshold" . 
The strength of electron-phonon coupling in graphene 
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FIG. 1. (left) The electron energy spectrum in graphene, 
exposed to a perpendicular quantizing magnetic field. The 
dashed line represents the phonon emission threshold energy 
e^ n . Here it corresponds to the Landau level e M „ with /i = — 1 
and n — 1. In contrast to the prediction from Ref. 1201 no 
state corresponds exactly to the energy e^ n . Instead, the 
true spectrum includes a sequence of electron-phonon bound 
states, which coagulates to the threshold e^ n . The most dis- 
tant bound states with the binding energies W\,2 of opposite 
signs are shown in a larger scale by the bold solid lines below 
and above elu- (right) Schematic diagram of the absorption 
spectrum in graphene. The cyclotron-phonon resonance is 
governed by the peaks corresponding to the electron-phonon 
bound states, which constitute an asymmetric doublet around 
v = UJ - 

a ~ 0.02 [23 while the phonon energy ojq ps 196 meV 
and we find that the binding energy should be of the or- 
der of Wgr = 4 meV. Because of its different functional 
dependence on the coupling constant a, the the bind- 
ing energy in graphene is much larger than that in bulk 



2 



samples [23] or in graphene samples in the absence of 
a magnetic field |25j . In conventional two dimensional 
electron systems the binding energy scale in a magnetic 
field has the same functional dependence on a [55] , how- 
ever due to the larger phonon energy, W gr is substan- 
tially larger in graphene than binding energies in InSb 
(W 2D = 0.5 meV) and GaAs (W 2D = 2.5 meV) struc- 
tures, extensively studied in experiment (27j [28]. Our 
calculations show that the dependence of the binding en- 
ergy in graphene on the magnetic field is weaker than that 
in conventional systems and for 0b ~ 0o actual binding 
energies vary within 5—15 meV that is fully measurable. 
Moreover in graphene in such moderate fields, the bound 
states can be located between the zero and first Landau 
levels with fi = +1 (see Fig. [TJleft)) and this creates 
an additional possibility that the splitting of the absorp- 
tion peak can also be observed in combined resonance 
absorption with a participation of the zero Landau level. 
In this case the light frequency can substantially differ 
from the frequency of transverse optical phonons, near 
which the lattice reflection is important [551 130] ■ This 
in its turn facilitates the experimental observation of the 
bound states in graphene. 

In experiment the bound states will manifest them- 
selves in the absorption spectrum of the cyclotron- 
phonon resonance |31l 132] , Upon absorbing a photon 
of frequency close to v — e^ n ' — e^ m + ujq, an electron 
is transferred from the initial Landau level \in to the fi- 
nal state i/n' and simultaneously it creates an optical 
phonon, which is bound to the electron in its final state. 
This process results in a fine structure of the cyclotron- 
phonon resonance spectrum S(v), which for n — n' and 
fx = fi' is shown in Fig. [IJ right). Instead of a single 
delta- function peak at v = ljq, the absorption spectrum 
constitutes an asymmetric doublet around the phonon 
frequency uiq with a splitting of the peak determined by 

V-LOq- Wgr [33]. 

Threshold approximation and diagrammatic technique 
We calculate the spectrum of the electron-phonon 
bound states from the poles of the electron-phonon scat- 
tering amplitude S(e) with respect to the energy pa- 
rameter e, corresponding to the total energy of the elec- 
tron and phonon [33]. Despite the weak bare electron- 
phonon coupling, the effective electron-phonon interac- 
tion is strong in the energy range close to the phonon 
emission threshold e^„ where the density of final states 
of the electron+phonon system diverges. Therefore, the 
calculation of by means of perturbation theory is not 
applicable [see the review papers in Ref . [3D] . The exact 
electron-phonon scattering amplitude satisfies the Dyson 
type integral equation, drawn in Fig. [2] in which we have 
explicitly separated a "dangerous" intersection near the 
threshold with respect to one electron and one phonon 
line. Such an intersection corresponds to the emission of 
an almost real optical phonon and results in a singular 
term, which is responsible for the enhancement of the ef- 
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FIG. 2. Equations for the electron-phonon scattering am- 
plitude, (top) The exact scattering amplitude E(e) with a 
" dangerous" intersection with respect to one electron and one 
phonon line, (bottom) The irreducible four vertex part D(e), 
approximated by the bare amplitude. The solid and dashed 
lines correspond to the electron and phonon Green functions, 
the bold dots to the bare electron-phonon vertices. The ver- 
tical dotted line shows the "dangerous" intersection. 

fective electron-phonon coupling and for the formation of 
the electron-phonon bound states. 

The diagrammatic equations in Fig. [2] are depicted in 
a gauge invariant diagrammatic technique (details about 
this technique can be found in the Supplementary infor- 
mation and in Refs. [30] and I35|) . The filled and open 
squares correspond to the exact, £ Mrii/ y„/ (e|q, q'), and 
the bare, □ /in!AI '„'(e|q, q'), electron-phonon scattering 
amplitudes, the internal solid and dashed lines to the 
exact electron and phonon Green functions, G At „(e) and 
D(p, q). All these quantities do not depend on the elec- 
tron gauge non-invariant quantum number k y , i.e. on 
the choice of the vector potential gauge. In the "dan- 
gerous" intersection in the second term of the rhs of the 
diagrammatic equation for E(e) in Fig. [2j we can at low 
temperatures replace the exact phonon Green function 
by a free phonon propagator. After taking the integra- 
tion over the phonon energetic parameter w, it is replaced 
by the phonon energy u/q. Then, the integral equation for 
E(e) in Fig.[2]can be written analytically as 

S M „ )(tl '„'(e|q,q') = □ Am , Al '„'(q,q') + D^ n<m (q_, C0.) 

xG m (s - uj q ) exp (-iAx) ^p.n,^' n ' (e|q, q') • 

Here the exponential factor appears due to the gauge 
non-invariant part of the electron-phonon vertex (see 
Eq. [3] and the Supplementary information for details). 
The phase Ax = ^[q — q, q — q']/2 where the square 
brackets denote the out of plane component of the vector 
product. 

In the energy range of our interest e ~ £^„ the energy 
parameter e — uj of the electron Green function in the 
second term of the rhs of Eq. ([T]) is far from the threshold 
e[\ni therefore the exact electron Green function can be 
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taken without allowance for interaction with phonons and 
be replaced by the bare function G^ n (e — u>o)- I n the 
Landau representation G°„(e) is diagonal and is given 
by its usual pole form 



G> n (e) = (e - e^ n + iO ■ sgas)' 



(2) 



Therefore, in the sum in Eq. over the electron internal 
quantum numbers /2, n only the singular term with ft = fi 
and n — n should be retained in the neighborhood of the 
phonon emission threshold e ~ e c . 

The central quantity in Eq. (fTJ) is the bare amplitude 
n^n,^fi(q, q), which has no "dangerous" intersection with 
respect to one phonon and one electron line and there- 
fore can be calculated within perturbation theory. In 
the lowest order in a, the two diagrams shown in Fig. [2] 
contribute to □ A1Ji ,^ fi (q, q). The bold dots represent the 
bare electron-phonon vertices 7(q), which we calculate 
in Supplementary information. For longitudinal optical 
phonons we find 

l^nk y ^nik' y {q) = \fa-jr-$v y ,k y +q x exp [i^ q (n - n')} 

x exp [iq x (k y + k' y )/2] P^ n)/i / n ' (i) (3) 

where C is the normalization length, (f> — 4>q the polar 
angles of the vector q. The gauge-invariant part of the 
vertex is 



Pfin^u'n' (0 ( 1) 



[//< 



'-l(t) ~ M<2n-l,n' (*)] 
(4) 

with Co = 1 and c„ = \j\[2 for n > 1. For n > m 
the form factor Q nro (t) = V m! /n!e- t/2 i (n - m)/2 ££r m (i) 
where L^(t) is the Laguerre polynomial and t = q 2 l 2 B /2. 
For rt < m, one can use the identity Q n m(t) = 
(—l) n ~ m Qmn(t)- For either n = or m = the form 
factor Q nm (t) with a negative index is identically zero, 
which means that in graphene the electron-phonon ma- 
trix elements between the zero-energy Landau states van- 
ish identically, P^o^'oC*) = 0. 

As seen from Eq. ([3|, the gauge non- invariant part of 
the electron-phonon interaction vertex in graphene is the 



same as the one in conventional electron systems with 
parabolic dispersion [36]. Consequently, it allows us to 
apply the rules of the gauge invariant diagrammatic tech- 
nique for conventional systems from Ref. [55] when calcu- 
lating the electron-phonon scattering amplitude in 
graphene. We just replace the gauge invariant form fac- 
tor Q nm (t) by P^n.n'n'it)- 

Dispersion equation of electron-phonon bound states - 
To find the spectrum of the bound states it suffices to 
consider the amplitude S(e) with the fi — [i! and n = n' 
external electron lines. After the above simplifications 
we can also eliminate the phase factor Ax in Eq. |TJ) 
by introducing a symmetric definition for the scatter- 
ing amplitude E'(e|q, q') = exp ^— i^f- [q, q']J S(e|q, q') 
(similarly for D'(q, q')). These primed quantities will 
depend only on the difference of the polar angles 
4> = <f>q — <pq' of the vectors q and q'. There- 
fore, we find that the Fourier components S^ n (£;tt') = 

/ Q 27r dcj>/ '(2tt) exp(— ^0)S^ m (e|q, q') of the new amplitudes 
with different I = 0, ±1, ±2, . . . satisfy the independent 
integral equations 



RUs\t,t') = KUt,t') 



+ A M „ (e) / diKi n (t, T)Rln (#> 



(5) 



Here we define the dimensionless amplitudes, 
<n(e|M') = (Tr~f/au )^ n (e\t,t') and i^ n (e|M') = 
{■K-i/au )U'l n {t,t') where 7 = g s g v C 2 /2tt£ 2 is the 
capacitance of the Landau level allowing for two spin 
orientations, g s = 2, and for the valley degeneracy, 
g v = 2. We introduce also the function 



A M „(e) 



CJ 



4-7T e — e c . 



(6) 



which determines the effective electron-phonon interac- 
tion in the energy region near the threshold e^ n with the 
renormalized coupling strength a = awg/wo- As seen 
the amplitude R l (e\t, t') is the resolvent of the inho- 
mogenous Fredholm integral equation ^ with the kernel 
K l (t, t'), which we can write explicitly as 



E (- 



iy- n p, 



fin, ism \ 



J l-\-m— n (2> 



$L . m. — r 



(yyjm, — \i\]n) + cr (y^Jrri — li^fn) — o 



(7) 



Here a = ujq/ujb, Ji(t) is the Bessel function, and Si^ m 
is the Kronecker delta. Hence the resolvent R l (e\t, t') 
has a pole in e when A^ n (e) = A' r where \\ inr are 
eigenvalues of the kernel Q , numbered by the index r — 
1,2,... for a given I. Thus, the energies of the electron- 



phonon bound states are given by 



c fxnr c fin \[ — c un v " unr V°7 

where W l ^ nr defines the binding energy of the bound 
states, which refer to the threshold e^ n . As far as there 
is no continuum in the bare spectrum, we find that the 
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new branches of the spectrum are not damped and ap- 
pear in the spectrum both below (Wh nr > 0) and above 
(W,Lr < °) the "threshold" energy E c „ n , i.e. the binding 
energies are determined by the kernel eigenvalues A' r 
of both signs. The bound states are identified by the to- 

TABLE I. Eigenvalues of K l r = (—1)' k,7 ; in strong magnetic 



fields. 
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0.97 


-5.08 
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-23.49 


53.68 


-126.5 
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-3.27 


5.19 


-10.24 


21.22 


-46.19 


103.4 


2 


3.42 


-6.38 


13.23 


-28.68 


64.14 


-146.6 



tal angular momentum I of the electron rotation around 
the magnetic field, which takes both non-negative and 
negative integer values. For a given / there exists a se- 
quence of bound states, labeled by the index r. The 
latter emerges due to the integral nature of the equation 
for the scattering amplitude E(e). Because the disper- 
sion of the optical phonon is negligibly small, the energy 
in the "dangerous" intersection does not depend on the 
phonon momentum q: all intersections with different q 
are singular and phonons with any q participate in the 
formation of the bound states. 

Binding energies of electron-phonon bound states - 
The analysis of the eigenvalue problem for the kernel 
([7]) shows that in high magnetic fields (a <C 1) the lead- 
ing order in tr _1 contribution to the sum in ([7| comes 
from the term v = [i and m = n. Hence, in this regime 
the effective electron-phonon coupling for electrons from 
the n = 1 Landau level, P^x tlJ ,i{t), (recall for n = 
the electron-phonon vertex P^,o^o(t) = 0) determines the 
kernel 

KUt,t') = -Qi (t)Q 10 (t') U(2Vtt 7 ) - <5 ; , 



(9) 



which is independent of the chirality /1. The eigen- 
values even of this simplified kernel cannot be found 
analytically. However, the parameter a enters in the 
kernel as a factor so its eigenvalues are given by 
X l ±lr = aK l r where K r are the eigenvalues of the ker- 
nel Vtt'exp [-{t + f')/2] J ^ (2^/tt 7 ) - 8 h0 , which we cal- 
culate numerically and summarize in Table [T] Thus, in 
this strong field regime the binding energies are given by 
Wj. = au>o I (a 2 K l r ) , which increase linearly with B and 
decrease rapidly with r, asymmetrically on both sides 
of the threshold. In weak magnetic fields (cr 3> 1) the 
spacing (y/m + 1 — y/m)u)B ~ coo/2a 2 between adjacent 
Landau levels m and m + 1 around the phonon emission 
threshold becomes smaller than the scale of the bind- 
ing energy W ~ auiQ. Therefore, in graphene already 
for l/y/2a < a ~ ^/m Landau levels form a quasi- 
contimmm, which hinders the formation of the bound 
states in such weak fields. 

For arbitrary magnetic fields we calculate the binding 
energy of bound states numerically from the full kernel 



30 
£ 25 

4 20 

7 10 
5 



B 2 



1 











1 ' 

1 

1 


1 

/ / 
/ / 

/ / 








/ / 








1 / / 
















J 













1.0 1.5 2.0 3.0 5.0 
B[B ] 



FIG. 3. The binding energies W® nr for the bound states, 
which refer to the threshold e^ n with n = 1 and fi = — 1. The 
solid and dashed curves correspond to the r — 1 and 2 bound 
states, which are most distant from the threshold. 



([7]). Near the resonances e^ n ' — e^ n = ±oj between 
the bare graphene Landau levels, electron-phonon hy- 
brid states are formed in the spectrum, which have a 
much larger binding energy scale of the order of ^/oluiq. 
These resonance states can be obtained within perturba- 
tion theory from the poles of the single- or two-particle 
Green function (but from the poles of the scattering 
amplitude) and determine the structure of the magnc- 
tophonon resonance. This effect has been recently stud- 
ied in Refs. l9l and H0l and is here excluded from the con- 
sideration. Far from the resonances the spectrum is com- 
pletely governed by the electron-phonon bound states 
([8]), which determine the fine structure of the cyclotron- 
phonon resonance (c/. Fig. [l] (right)). In Fig. [3] we plot 
the binding energies W_x\ T as a function of the mag- 
netic field strength B for the most distant r = 1 and 2 
electron-phonon bound states with / = 0, referring to the 
threshold e c _h. The divergences of the binding energy 
correspond to the hybrid states at the resonance fields 
B = B , B = Bi, and B = B 2 . Here B Q 22.04 T, 
Bi fa 1.87-Bo and B 2 « 5.83£?o are given, respectively 
by the resonance conditions ujq = ojb, ^0 = (v3 — 
and ojq = (y/2 — It is seen that for intermediate 

magnetic fields far from the resonances Bq 7 Bi, and B 2 , 
the actual binding energies of the electron-phonon bound 
states vary approximate from 5 to 15 meV and are mea- 
surable in light absorption experiments. 

In conclusion, we have calculated the electronic spec- 
trum in a monolayer graphene in a perpendicular mag- 
netic field. It is shown that the true spectrum near 
the optical phonon emission threshold consists of new 
branches of electron-phonon bound states, which will 
manifest themselves in the absorption spectrum as an 
asymmetric doublet around the threshold energy. 
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